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Abstract 

We study the gravitational effect of non-self-annihilating dark matter on compact stellar objects. 
The self-interaction of condensate dark matter can give high accretion rate of dark matter onto 
stars. Phase transition to condensation state takes place when the dark matter density exceeds 
the critical value. A compact degenerate dark matter core is developed and alter the structure and 
stability of the stellar objects. Condensate dark matter admixed neutron stars is studied through 
the two-fluid TOV equation. The existence of condensate dark matter deforms the mass-radius 
relation of neutron stars and lower their maximum baryonic masses and radii. The possible effects 
on the Gamma-ray Burst rate in high redshift are discussed. 
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I. INTRODUCTION 



In the ACDM cosmology model, dark matter contributes 24% of total energy of the 
universe. Dark matter is regarded as the majority of matter in the universe, and it has sig- 
nificance effect on large-scale structure formation and evolution. Astrophysical observations 
provided indirect means to probe dark matter and study its effects. 

Dark matter may also have effects on individual stellar objects. Steigman et al. [l| 
put forward the idea that capture of WIMP particle would affect the stellar structure and 

n n 

evolution. Spergel and Press |2j and Faulkner and Gilliland [3[ employed the idea for an 
explanation of the solar neutrino problem. Recently, there are growing interests in this 
field. Self-annihilating dark matter provide extra energy for stellar objects, which have been 



extensively studied by 



by 



studied by jg], jo], [10|. The impact on the evolution path of main sequence star is studiec 



1 1 1 



7]. In particular its impact on first-generation stars is 



fl, 



111 ]. 12J. For non self-annihilating dark matter, its impact is studied by [13] 
for main sequence stars and for neutron stars by 16|], 17], 18( 19| for different dark matter 
models. 

In this paper, we study the gravitational effect of non-self-annihilating condensate dark 
matter on stellar objects. The Bose-Einstein condensation may come from the bosonic 
features of dark matter models, e.g . axions, Nambu-Goldston bosons or supersymmetric 
partner of fermions suggested by J20 [| . 



21 |.|22| . Some previous studies of boson star structure 



and dynamics can be found in 23[, 24], 25|,26[. Condensate dark matter has an equation of 



state P v 



2irh 2 L 



p: 2 , where l a is the scattering length and m x is the dark matter particle mass. 



Condensate dark matter is able to form gravitationally bound objects after the condensation 
takes place. The properties, structure and stability of such objects is studied by {27]. Phase 
transition to condensation can occur under either one of the following two conditions: when 
the temperature cools below critical value or when the density exceeds the critical value j^ . 
In this paper we consider the following picture: the deep gravitational well of a star accretes 
dark matter, when phase transition takes place, dark matter forms a compact degenerate 
dark matter core inside the star and may causee the star to collapse to compact object, 
i.e. white dwarf, neutron star or black hole. We examine conditions of the existence of 
such degenerate core can make the white dwarf or neutron star stable. If no stable neutron 
star is possible, the star will eventually collapse to a black hole. The accretion rate and 
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criteria for phase transition are presented in this paper. Following that, analytical and 
numerical study for the structure and stability of various stellar objects with condensate 
dark matter component are presented. The parameters for condensate dark matter are 
l a = lfm, m x = lGeV as fiducial values. 



II. ACCRETION OF DARK MATTER ONTO STARS 



The self-interaction of condensate dark matter enables the dark matter particles to trans- 
fer energy to each other. Through interacting and losing energy to dark matter particles 
already inside the star, incoming dark matter particles will be captured by the star and 
stay inside it. The calculation of accretion rate of dark matter particle follows from Gould 
while self-interaction is considered the interaction between dark matter and baryons is 



ignored. The self-interaction cross section of dark matter is determined by 

a xx = Anil (1) 

Assuming spherically symmetric gravitational field for stars, the escaping velocity is v. 
Dark matter has velocity distribution f(u)du. Q~(w) is defined as the rate per unit time 
that a dark matter particle with velocity w will be scattered to a velocity less than v when 
collided with materials. The flux of dark matter particles that goes inward is 

^f(u)ududcos 2 9, (2) 

where 9 is the angle relative to the radial direction, < 9 < |. The probability for a dark 
matter particle to be scattered to velocity less than v is 

n«(«0^, (3) 



where 

dl 2 



i^ 2 

rw 



'drQ(rw-J). (4) 

w w l \rw/ J 

J is the angular momentum of an incoming dark matter particle. Q(rw — J) is the step 
function. After integration over whole angular momentum space, the total capture rate per 
unit shell volume is 

dC f°° , f(u) „ , . 
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Let the dark matter particles with mass m x and number density n m in the star. From 
kinematics, the energy loss in a collision, ^ lies in the interval 

AE 

< -g- < 1, (6) 

The net capture of a dark matter particle requires 

w 2 — v 2 u 2 AE v 2 

= — < < — . (7) 

w 2 w 2 E w 2 

The total probability of capturing a dark matter particle is 

= a xx n m w—^Q (v 2 - u 2 ) . (8) 

Taking the incoming dark matter particle follow the Maxwell-Boltzmann distribution with 
number density n ont , 

f(u)du = n out ^=x 2 exp(-x 2 )dx, (9) 
\/7T 



where x is the dimensionless velocity defined by 



2 m x u 2 3u 2 
x 2 = —A — = — (10) 
2k B T x 2v 2 v ; 



with v being the average velocity of dark matter particles in galaxies. 
The capture rate per unit shell volume is 



dC /fi\^ v 2 



where A — Generally A ^> 1, the mass accretion rate is given by 



M x = j m x dC = ^M x (t)a xx n°^. (12) 
So the total mass accreted at time t is 

M x (t) = M exp(^a xx n OUt ^)- (13) 

The accreted dark matter mass grows exponentially. However, it is only valid when the 

probability term ([3]) is smaller than 1 and shall never exceeds 1. For a star with radius i?*, 
an estimation of when the probability reaches 1 can be given by 

n in a xx R* = 1. (14) 



When the probability term ([3]) equals 1, all incoming dark matter is accreted onto the star. 
So the capture rate is simply integration of fl2J) over the star surface. 



/oo 2,71 
m x f(u)udu = — 



2 / — \ / out 



2 X W ' WM ^ ' Tofc) (250^) ( l<w£ V /cm» ' ' (15) 



Hence, the total accretion rate is given by 



M x 



l M x a xx n^- n™a xx R* < 1, 
'f Rln out vm x n in a xx R* > 1. 

In our calculation, dark matter mass will first grow exponentially. After reaching a certain 
mass, the accretion rate becomes constant. 

The total dark matter mass accreted is given by 

2 



M x (t) 



M x0 exp 



c QO ( la \ 2 ( n°"A ( V \ ( y \ ( 
°- 00 Klim) yW® J ^516km/sy ^250km/sy ^ 



t 

lyr 



n m a xx R* < 1, 



M xl + 2 x 1O- 1O M (j^) 2 (^) (gS) (^) n'V xx ^ > 1. 



M x q is the initial mass of dark matter inside the star before accretion and M x i is the 
mass of dark matter inside the star when the probability term ([3]) reaches 1. With typical 
values of parameters of condensate dark matter, the constant accretion rate is reached very 
quickly. For a star with solar mass and radius 10 11 cm, let 77° ut = 10 10 , = m x n ont , the 
probability takes place when n m w 10 14 . So the exponential accretion lasts for lo |gg w 1.5yr. 
The lasting time of exponential accretion is inverse proportional to the logarithm of initial 
mass, which will be very small compared to the stellar lifetime. Therefore the constant 
accretion stage actually determines that the amount of dark matter can be accreted within 
the stellar lifetime. The total dark matter accreted within the stellar life is proportional to 
the ambient dark matter density, if the intrinsic properties of dark matter are fixed. 

The density of incoming dark matter particles near the star constrains the accretion rate. 
In the early universe the ambient dark matter density is much higher than today. When the 
early stars are formed, the creation of deep gravitational well in the dark matter halo makes 
it to further contract and increase the density of dark matter near stars to be even higher. 

The increase of dark matter density in star formation can be calculated by using the 
adiabatic approximation. The adiabatic invariant is rM(r) where M(r) denotes the total 



5 



mass enclosed within radius r. This leads to the adiabatic equation given by Blumenthal et 



al. 



29] 



M (r ) 

M(n) ri = = (M x ( ri ) + M b (r f ))r f . 



(16) 

1 _ Jb 

M(r), M x (r) and M^(r) are respectively the total mass, dark matter mass and baryon mass 
inside radius r. ri,rf denote the radius before and after the contraction inside which dark 
matter of total mass M x (rj) is enclosed, fa = 0.15 is the mass fraction of baryon in the 
initial dark matter halo. 

Fig |T] il lustrates the contracted dark matter profile for different initial profiles: NFW 
profile 30|, isothermal sphere 3l| and Burkert profile 32|. Dots are from the simulation 
results by Abel et al. 33[ The parameters of each initial profile are tuned to fit the simulation 
results. The best fit line is from the NFW profile. The parameter for the dark matter halo 
is chosen to have total mass 7 x 1O 5 M and the concentration parameter as in 30| c = 2. 



This result agrees with the calculation done by Freese et al. 



34] 




FIG. 1: Comparison Between Adiabatic Contracted Density Profile of Different Initial Dark Matter 
Distribution and Simulation Result by Abel et al. 33j]: Dots are simulation results. 



The adiabatic contracted dark matter density profile fits the simulation result quite well 
from O.Olpc outward, but how close to the core this profile is still valid is unclear, as near 
the core the strong diffusive process in star formation may break the adiabatic condition of 
dark matter. 
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III. PHASE TRANSITION OF DARK MATTER INSIDE THE STAR 



Condensate dark matter particles will undergo Bose-Einstein condensation inside the star 
when its density is above the critical one. The critical density for this kind of phase transition 
is given by [27] 



- 7.327 xlO'f f',,, ) f^m" g /W, (17) 



27rh 2 l a \9 x 10 14 cm 2 /s 2 J VI GeV/ V 1 fm , 

where a 2 v is the velocity dispersion of normal dark matter. 

The distribution of dark matter particle inside the star is beyond current knowledge, 
while it might be possible to extrapolate some from the distribution without the gravitational 
potential. Simulations and self-similar solutions give that dark matter density of dark matter 
profile will follow the power-law distribution p oc r~ a (cf [30|, |35[, |36j, (37J). The exact 
value of a is still unclear, but most studies show the value lies in the range 1 < a < 3. With 
the existence of extra gravitational well due to the baryonic star, the dark matter profile 
will have a steeper profile, for more dark matter will be dragged into the core by the extra 
gravity. 

In a very small central region that dark matter density may exceed the critical density, 
the dark matter within this region will become condensate but they can only form some 
microscopic droplets before they reach the minimum mass of the gravitationally bounded 
object (cf (27)]), whose exact size depends on the detail of the exact interaction of dark 
matter and their surface energy. When dark matter continue to accrete into the core of the 
star, the number of microscopic droplets increase. When the total mass of these droplets 
is sufficiently large, they become the gravitationally-bound object. The minimal mass of 
stable condensate dark matter object is obtained from j^] by setting the central density 
to be the lowest one which is the critical density of phase transition. With our selection of 
dark matter parameters, the minimal mass is M Ximin = 2.296 x iO~ 5 M . 

Taking the total dark matter mass inside the star to be M x (t) = J M x dt, to be the 
stellar mass, and the dark matter density profile to follow 

Px = Pxfi (jj-) > ( 18 ) 
the dark matter mass in the region where density is above critical value is given by 

3-a 

Pxfi X 



^X,cen = M x [f±) . (19) 



The criteria for the condensate dark matter core to form is 



-^Xi cen ^ -^Xi m i n ' 



(20) 



IV. FORMULATION FOR STARS WITH DARK MATTER COMPONENTS 



We assume that dark matter only interact with baryons through gravity and study the 
static structure of such stars with dark matter components. We adopt the two-fluid formu- 



lation by 



16], which is first introduced by [17J. This formulation is a separation of TOV 



equation motivated by the similarity of structure equations between the relativistic and 
Newtonian ones. There are two equations for the balance between gravity and pressure for 
baryons and dark matter separately. 



,dP*(r) 



dr 



,dP, 



DM in 



dr 



-GM(r) 
-GM(r) 



1 + 
1 + 



47rr 3 P(r) 
M(r)c 2 

47rr 3 P(r) 
M(r)c 2 



1 - 



2GM(r) 

c 2 r 
2GM(r) 

c 2 r 



-i 



Pb{t) 



Pvmir) 



(21) 



where P(r) = Pe(r) + PdmM and M(r) = M B (r) + M DM (r). With two more equations of 
mass continuity of baryon and dark matter separately, 

dM B (r) 



dr 

dM DM (r) 
dr 



4:irr 2 p B (r), 
4nr 2 p DM (r) 



(22) 



the four equations give a complete set of equations for the structure of stars with dark matter 
component. 

Even though there are two-fluid formulation available from first principle calculation 



they give identical numerical results for static stellar cases 39J. Therefore, we take this simple 
approach to study the static stellar structure with dark matter component. 

Once dark matter condensation takes place, the gravitationally bound object formed has 
density at least 7.327 x 10 9 g/cm 3 27J, much greater than the onset of electron degeneracy 
10 4 — 10 5 g/cm 3 ) 40J and central density for most white dwarfs (10 4 g/cm 3 — 10 10 g/cm 3 ) 
4=lj | . For normal stars or white dwarfs, the baryonic component inside the dark matter core 
is negligible. The normal dark matter interacting with the condensate core will also rest 
to the ground state, become condensate and accrete onto the core. Therefore, we assume 
that all dark matter inside the star are concentrated in the degenerate core. Under this 



S 



approximation, the star has a pure dark matter core and outside the core is pure baryonic 
matter. We treat the pure dark matter core as initial conditions for the normal TOV 
equations and integrate the equations for baryonic matter from the surface of the dark 
matter core. 

The normal TOV equations for stars is given by 
dM 



dr 

dP 

dr 



Airr p(r) 



G 



M{r) 



4irr 3 P 



2GM(r) 
c 2 r 



-i 



(23) 



The dark matter core has radius r c and mass M c which can be determined from [27J and the 
density of normal baryonic matter at the surface of the core is po- The initial conditions for 
the TOV equations for the star with dark matter core are when r = r c , M = M c , p = po- 
Using the parametrization below 

P P 







m 



m. 



a 



Cx c 



Mc—y/AnGpo, £ c = -y/AnGpo, 

Cx C 



the TOV equations reduce to the following form 42] 

-Ian 



dm 
d£ 
da 
d£ 

0(£ c ) 



m) 



e(i - f ) 



1, m(£ c 



m c . 



(24) 



(25) 



When £ — > £ c , the equation leads to 



do- (1 + cr)(m c + cr£; 



C 2 (l - f-) 



(26) 



In order to have physical solution, ^| < must be satisfied. This constraint leads to 



2—^ < 1, 



equivalently 



2GM r 



c*r c 



< I. 



(27) 



(28) 
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This constraint means that the core cannot directly collapse to a black hole. 

Another constraint comes from the non-relativistic limit of Eq.f l25|) in the region where the 
gravity of dark matter core dominates. So the mass and self-gravitation of normal baryonic 
matter are neglected, m « m c . In the non-relativistic limit, er <C 1, the TOV equations 
become 

" > = < 29 > 

Integrating both sides, since 9 < 1 the left hand side gives 

/ ~da>a(Q-a(C). (30) 
The right hand side can be integrated analytically, 



£ ui l 1 — 2mc 

,, c ^r^) de = 2 l0g T^- (31) 

Since the dark matter core is much smaller than the star, in regions £ 3> £ c , stable solution 
requires the pressure to be always positive 

^ c )>4 1 °g( 1 -^ £ )- ( 32 ) 

Since in Eq.( |30|) the integration takes 9 = 1, this analysis gives sufficient condition of p to 
have stable solutions. 

In case that baryonic matter has a polytropic equation of state P = Kp 1+1 ^ n , Eq. fl30p 
can be integrated analytically. In addition, if the core potential is non-relativistic 2 ^^ c 1, 
the constraint reduces to the form given by \a2\ 

> , GM \ ■ (33) 
c 1 (1 + n)c 2 r c 

This constraint gives the minimum value of po to have stable solutions. 



Po,, 



'1 + n)Kr c 



(34) 



V. APPLICATIONS 



Since normal stars have very low density, applying the analysis in Section [IV] to them 
gives that, once the condensate dark matter core is formed, the star is unstable except 
for extremely light dark matter core. The core of the star is expected to collapse to form 

10 



compact objects, i.e. white dwarfs or neutron stars. Instead of considering the actually 
collpase processes we only focus to study the stability of the static structure of white dwarfs 
and neutron stars with dark matter component. 



A. White Dwarfs 

As stated in Section IIVI the central density of most white dwarfs is much smaller than 
condensate dark matter and baryon component is negligible inside the dark matter core. 
We take the condensate dark matter core as initial conditions for the TOV equation and 
integrate the equation from the core surface. The equation of state (EOS) for baryonic 
component is taken to be the ideal fermi gas. This EOS is only valid below Neutron drip 



4 x 10 11 g/cm 3 40|) and neutronisation threshold (~ 1.22 x 10 7 g/cm 3 for inverse beta 



decay 40]). For high central density white dwarfs (> 10 7 g/cm 3 ), its equation of state has 



been calculated by [43[. When the central density becomes comparable to that of condensate 
dark matter, we need to solve the two-fluid equation. The integration of two-fluid equation 
is incorporated in the next section of neutron star where we use the BBPS EOS for low 
density region. In this section the upper limit of central density p is set to be 10 10 g/cm 3 
as fiducial value. 

Fig J2] illustrates the comparison of the relation between total mass (baryon mass plus dark 
matter mass) and central density for white dwarf with a core of O.O1M and that without 
a core. The existence of the core contributes extra gravity that causes the white dwarf to 
contract, resulting to a increase of central density and decrease of maximum mass. 

Figj3] compares the relation between total mass and radius for white dwarfs with various 
dark matter core mass. The increase of dark matter mass results in both smaller radius and 
smaller mass. For dark mass smaller than O.O4M , there is an instability of white dwarf 

dM 
dpo 



when ^p- < 0. This criteria gives unstable fundamental mode of radial oscillation 40 |. 



However, this elementary stability analysis shall be taken with care, since ^ > is not 
sufficient to conclude the white dwarf is stable. A proper stability analysis shall include 
all radial eigenmodes which is beyond the scope of this paper. But any white dwarf with 
higher mass will be definitely unstable. We determine maximum mass by this instability. 
When the core mass keeps growing, there is no instability of fundamental radial mode, the 
maximum mass is determined when the central density of white dwarf reaches maximum 
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FIG. 2: Relation Between Baryon Mass and po for White Dwarf with and without a Dark Matte 
Core of 0.01M Q 
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FIG. 3: Relation Between Total Mass and Radius for White Dwarf with Various Dark Matter Core 
Mass 

FigJH shows the relation between maximum baryon mass and core mass. The increase of 
core mass reduces maximum baryon mass. There is a sharp decrease when the core mass 
exceeds 0.04M Q . In this case the maximum baryon mass is not determined by stability 
consideration any more but the constraint on central density. When the core mass is even 
higher, over O.O6M , there is no stable solution for white dwarfs. In this case the white 
dwarf most likely collapses to neutron stars. We want to remark that although in our 
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calculation we assume Fermi gas EOS for the entire white dwarf, realistically the surface 
boundary of a white dwarf is defined when the density go to zero. Obviously when density 
is sufficiently low the mass is no longer degenerate therefore the Fermi gas assumption is 
not valid. However the mass in the region where the Fermi gas approximation breaks down 
is very small therefore it does not affect our analysis. 
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FIG. 4: Relation Between Maximum Baryon Mass of White Dwarf and Dark Matter Core Mass 



B. Neutron Stars 



In case of neutron stars, the mass of baryonic component inside the dark matter core is 
not negligible any more. We have to solve the two-fluid equations, it should be noticed that 
in some cases, degenerate dark matter radius can be larger than the baryonic mass radius. 
We will illustrate this case in Fig. 6. 

We consider several popular types of EOS for the baryonic component, including UVU 



44], APR |45j, SLy |46|, FPS [47j, RMF-soft and RMF-stiff [48 



491 ] . For the low density part of EOS, we adopt the BBPS 43] . Fig. |5] shows the Mass 



Bombacil and Bombaci2 



Radius relation for these EOS's together with pure dark matter. 

The radius of such stellar object is determined when density becomes zero. The solu- 
tion of two-fluid equations of generally gives different radius for baryonic and dark matter 
components. Fig J6] gives two examples of stellar structure for baryonic and dark matter com- 
ponents using UVU EOS. We can see that dark matter radius can be much larger than that 
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FIG. 5: Mass-Radius Relation for Different EOS and Condensate Dark Matter 

of baryonic matter. As we assume dark matter does not interact with normal matter except 
through gravity, we take the radius when the baryonic density becomes zero as the radius for 
the two-fluid neutron star from the observational consideration, while the observable mass 
is the total mass. This is very important. 

FigJT] illustrates the relation between total mass and baryonic mass radius with the dark 
matter mass fixed. When the dark matter mass is small, the increase of dark matter mass 
deforms the original mass-radius relation, reducing the mass, radius and maximum mass. 
When dark matter mass becomes dominant, the maximum total mass increases again. For 
some EOS (RMF-soft and Bombacil), the maximum mass can be higher than that of pure 
baryonic neutron star. The extreme case of a dark matter rich two-fluid star is that it 
becomes pure condensate dark matter. Condensate dark matter star has a maximum mass 
~ 1.8M , so the maximum mass will finally reach this value. For dark matter rich two- fluid 
stars, the maximum total mass increases with dark matter mass and in some cases exceeds 
the maximum mass of pure baryonic neutron star. 

FigJH] shows the relation between total mass and radius with dark matter energy ratio e 
fixed. The energy ratio e is defined by e = j-j—^m~ ■ The increase of e reduces the mass 
and radius while generally keeping the shape. When dark matter dominates, the shape 
transforms that of pure dark matter, but the radius is smaller. 

Figj9] shows the relation between maximum mass of baryonic component and maximum 
total mass with the dark matter mass. The increase of dark matter mass always reduces the 
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FIG. 6: Examples of Structure for Two-fluid Neutron Stars 

maximum baryonic mass, but the maximum total mass increases after dark matter becomes 
dominant. 

VI. DISCUSSION 



The dark matter density outside the star is a key factor for accretion. After recombina- 
tion, the dark matter density evolves proportional to (1 + z) 3 . In the early university the 
ambient density of dark matter is much higher and they can form high-density dark matter 
halos. Simulations reveal "cuspy" density profile of dark matter halos 30j. The density 



increases exponentially towards the halo center. Even though strong density cusps are not 
observationally favorable (cf. [501]). the existence of deep gravitational well (massive stars 
or black holes) will pull the dark matter and increase the density around them. We have 
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FIG. 7: Mass-radius Relation for Different EOS With Various Dark Matter Masses 



calculated the increase of density using Blumenthal's method 



motion of dark matter particles. Freese et al. 34| have calculated the change of density 



29j | assuming only circular 



using a general treatment of adiabatic compression by Young 



5l| considering the noncircu- 



1(3 



lar motion. The difference between the two methods are within a factor of 2 [34J. Though 
adiabatic contraction formally requires the particles' orbital time to be shorter than the 
collapse time, in practice the method works well beyond this limit The contracted dark 
matter density can be ~ 10 10 GeV/cm 3 at 10 -4 pc from the gravitational center 34j. We 
hope that more knowledge shall be obtained by future simulations and observations of small 
scale structures. Furthermore although the average dark matter density in solar system is 
very low, it is possible that much higher density dark matter may exist near the galactic 
center, e.g. 
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531 ]. so stars in this region may accrete enough dark matter. 



By assuming the dark matter density is sufficiently high around the stars,e.g. in the dark 
matter halos of the early universe or at the center of the galaxies, we study the gravitational 
effect of condensate dark matter on various stellar objects. The deep gravitational well 
created by the stars pulls dark matter onto it. The self- interaction of dark matter particles 
enables them to lose energy and rest inside the star. When the dark matter density exceeds 
the critical value, the phase transition to condensation state takes place. Condensate dark 
matter develops a compact core and alter the structure and stability of the stellar objects. 
Normal stars will become unstable except for extremely light dark matter core. For white 
dwarf, the maximum dark matter core is about 0.06M Q . Condensate dark matter admixed 
neutron stars is studied through the two-fluid TOV equation. The existence of condensate 
dark matter has significant impact on the mass-radius relation and lower the maximum 
mass compared to neutron stars. Condensate dark matter admixed neutron star can result 
in compact objects with very small coordinate radius of visible object, which is determined 
by the baryonic matter. The existence of degenerate dark matter will also lower the mass 
threshold for stars to collapse to a black hole. This may lead to more GRBs caused by 
stellar collapse. Recent studies show that the rate of GRBs does not strictly follow the star 
brmation history but may be actually enhanced by some other mechanisms at high-redshift 



54 



561 ] . The standard collapsar model indicates that stars with mass larger than 30M Q 
can produce GRBs [57]. The existence of degenerate dark matter core can reduce the mass 
threshold of stars that can produce GRBs. If 3OM is the threshold for forming GRBs from 
stars, the threshold may be down to 2OM after including gravitational effect of condensate 



dark matter effect. For a Salpeter initial 
cutoffs mi = 0.1M Q and m u = 1OOM 



mass function of stars with lower and upper mass 



5N 



the percentage of stars that can die as black 



holes increases by 90%. So the star mass threshold of GRBs formation decreased due to 



17 



dark matter effect could be an alternative solution to the excess of the high-redshift GRB 
rate. 

Leung, Chu and Lin 



also studied the property of dark matter admixed neutron stars. 
We differ from them in several key respects. First, they assume dark matter as ideal Fermi 
gas but we study the Bose-Einstein condensate dark matter, which have very different EOSs. 
Second, Leung et al. 



did not consider how dark matter gets into the star, what are the 
key factors to determine the mass of dark matter inside the star and how dark matter is 
possible to rest inside the neutron star. We consider the accretion process through dark 



matter self-interaction from the surrounding halo. Finally, Leung et al. 
two-fluid model from first principle calculation. Though the numerica 



studied the 



16[ are identical 



results from their 
39[, their approach 



formulations and ours formulations given by [17J and 

can be used to study the time-dependent properties, like stellar oscillations 

We thank Dr. T. C. Harko and Prof. M.C Chu for useful discussion and suggestions and 
the anonymous referee for his very helpful suggestions and comments. 
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FIG. 9: Relation Between Maximum Mass of the Star and Dark Matter Mass. The upper panel 
is the maximum baryon mass vs dark matter and the lower panel is total mass (Baryon mass plus 
dark matter) vs dark matter. 
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